Critical temperature and superfluid density suppression in disordered high-T c cuprate 

superconductors 
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We argue that the standard Abrikosov-Gorkov (AG) type theory of T c in disordered d-wave super- 
conductors breaks down in short coherence length high-T c cuprates. Numerical calculations within 
the Bogoliubov-de Gennes formalism demonstrate that the correct description of such systems must 
allow for the spatial variation of the order parameter, which is strongly suppressed in the vicinity 
of impurities but mostly unaffected elsewhere. Suppression of T c is found to be significantly weaker 
than that predicted by the AG theory, in good agreement with experiment. 
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Sensitivity of the superfluid density and the critical 
temperature to moderate amounts of substitutional dis- 
order served as an early indicator of the unconventional 
nature of the order parameter in high-T c cuprate super- 
conductors. The theory of dirty d-wave superconduc- 
tors, pioneered in this context by Annett, Goldenfeld and 
Rcnn and by Hirschfeld and Goldenfeld [||, success- 
fully explained the observed T 2 dependence of the super- 
fluid density p s (T) in YBa 2 Cu 3 07_ < 5 (YBCO) films and 
single crystals, attributing it to the effect of the finite 
density of states at the Fermi level induced by impuri- 
ties, treated as unitary scatterers. Crossover to T-linear 
behavior was predicted for lower impurity content, which 
was later observed in very clean YBCO single crystals by 
Hardy et al. The theory was subsequently refined 
to provide a good quantitative description of the 
low-temperature behavior of p s (T). 

Similar models have been employed to predict the 
change of T c from its clean value T c0 due to disorder HQ , 
essentially by extending the standard Abrikosov-Gorkov 
(AG) theory M to the case of d-wave superconductors 
with scalar impurities. It has been noted that the ex- 
perimentally observed T c is much more robust than one 
would expect from these simple models, when measured 
against the corresponding change in p s (0). Figure [l] illus- 
trates this point by showing the experimental T c versus 
p s (0) for YBCO samples disordered by different types 
of disorder as obtained by various groups |P-]l3|| and the 
corresponding theoretical prediction Q . While there ex- 
ists a considerable spread in the experimental data, it 
is quite evident that the theory systematically overes- 
timates the suppression of T c , perhaps by as much as a 
factor of 2 in the cases of Refs. P, pd| , ^3[ . In order to rem- 
edy this situation, more realistic models have been con- 
sidered, taking into account strong coupling corrections 
within the Eliashberg formalism pj together with realis- 
tic band structures |1J] , proximity of the Fermi level to a 
van Hove point p| and the details of thepairing interac- 
tion within the spin- fluctuation model [M. While some 
improvement over the simple model can be achieved for 
carefully selected parameters, the discrepancy between 
theory and experiment remains in place, suggesting that 
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1. Normalized critical temperature versus the nor- 
malized zero temperature superfluid density. Experi- 
mental data were obtained by ^ mutual inductance J9| 
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methods traditionally employed to determine the sup- 
pression of T c by impurities are inadequate for the high- 
Tc cuprates. It has also been suggested that perhaps 
the effects of disorder on the transition from the rather 
peculiar and poorly understood normal state cannot be 
satisfactorily described within the framework of a simple 
BCS-like mean field theory @. 

In the present article we show that one can, in fact, for- 
mulate an adequate description within mean-field theory, 
if the spatial variations of the order parameter caused 
by random disorder are accounted for in a fully self- 
consistent manner. We present a simple argument for the 
breakdown of the conventional AG type theory (which 
enforces a uniform gap averaged over disorder) in su- 
perconductors with very short coherence lengths such as 
the high-T c cuprates. This argument is then given sub- 
stance by numerical calculations within the Bogoliubov- 
de Gennes (BdG) framework, carried out for a model d- 
wave superconductor with a random distribution of non- 
magnetic point impurities. Such calculations convinc- 
ingly demonstrate that the rapid drop in the predicted 
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T c is an artifact of AG theory resulting from spatial av- 
eraging of the order parameter. The true critical tem- 
perature obtained from a fully self-consistent solution of 
the BdG equations is always higher, in agreement with 
experiment. 

In a d-wave superconductor non-magnetic impurities 
are pair breaking and lead to suppression of T c . A simple 
modification of the original AG approach || yields an 
equation for T c of the familiar form [p|,f7||T4|] : 

K£)-*G + £0-*G)- (1 » 

where ip is a digamma function and a = l/2rT c o is a pair 
breaking parameter. In the limit of unitary scatterers, 
which is relevant for Zn and other solutes in YBCO 0, 
the scattering rate is given by l/2r = ni/ivN(0), with 
Hi being the number density of impurities and N(0) the 
normal-state density of states at the Fermi level. For a 
small we have 

T e AG 7T _ n, 

- 1 - - 1 - 4tv(o)T c0 ' (2) 

which is in fact an excellent approximation to the full 
solution of ([j]) for a < a c /3, where a c — 0.88191 is the 
critical pair breaking parameter \T c (a c ) = 0]. 

One of the crucial assumptions entering the derivation 
of Eq. (|l|) is that the position dependent order parameter 
A(r) can be replaced in the gap equation by its spatial 
average A. Such a procedure is valid when the gap varies 
over a length scale that is large compared to the average 
spacing between the impurities, U. In the vicinity of T c 
significant variations of the order parameter take place on 
the length scale set by the Ginzburg-Landau temperature 
dependent coherence length £ (T) , which is related to the 
low temperature BCS coherence length £o = Vf/nA by 
£(T) ~ 2/£ (l - T/Tc)- 1 / 2 , where v = 0.74 in a con- 
ventional s-wave superconductor with a spherical Fermi 
surface fll8| . For a d-wave superconductor v will be a 
different number of order unity, depending on the precise 
fc-space structure of the gap function. 

Consider the problem of a single unitary impurity in 
a d-wave superconductor. The order parameter will be 
strongly suppressed at the impurity site and it will re- 
cover its bulk value over the distance ~ £,(T) |l9| , ^0| . 
Since a single impurity cannot affect T c of a macroscopic 
sample, this coherence length is given by 

S(T)~ -T/Td,)- 1 / 2 . (3) 

For a small finite density of impurities with average spac- 
ing U ^> £o: except very close to T c , the areas of depressed 
order parameter around individual impurities will not 
overlap, and most of the sample will remain completely 
unaffected. Thus it is still reasonable to use Eq. (||) with 
the unperturbed T c0 for £(T). When these areas begin 
overlapping, i.e. when £(T) > k, the entire sample is af- 
fected in the sense that the order parameter is suppressed 



everywhere. The temperature at which this happens pro- 
vides a lower bound for the true T c of the sample. Using 
U = 2ao / ^TTrii (where ao is the ionic lattice spacing) this 
condition becomes 

^>l-j(^Vn, (4) 
T c o 4 \a Q J 

If T C AG predicted from Eq. (||) falls below this lower 
bound, one may conclude that the AG theory is not self- 
consistent, since the true coherence length at T = is 
smaller than the average distance between the impurities 
and the averaging over A(r) is not allowed. 

In order to get a rough idea of how restrictive this 
argument is, we note that the prefactor of rij in Eq. (|^) 
can be converted into a length scales ratio, 1/N(Q)T c q ~ 
(7r 2 /l-i3)(£o/a ) Q. Then, combining this with Eq. (g) 
one can derive a rough estimate for the range of validity 
of the AG theory of the form 

Vao > tt/1.13!/ 2 ~ 5. (5) 

Clearly, for conventional low-T c superconductors this 
condition is easily satisfied as £o is typically large com- 
pared to ao- For high-T c cuprates however, the situation 
is quite different since typically £o is several lattice spac- 
ings (£q w 4ao in YBCO). A more careful estimate using 
realistic values of T c q and N(0) for YBCO H], confirms 
that the result of our rough estimate (§) holds. This in- 
dicates that the usage of the AG theory for this material 
is probably not justified p2] |. 

While the argument presented above is admittedly 
crude, we believe that together with the experimental ev- 
idence discussed above it points to the necessity of study- 
ing the effects of spatial variations of the order param- 
eter in disordered short-coherence-length superconduc- 
tors. We now present results of a numerical calculation 
that strongly support the picture outlined above. 

We employ an extended Hubbard model on a square 
lattice with nearest neighbor attraction and on-site re- 
pulsion: 

H = - tJ2 4<7 C J<r ~ H ^2 + ^ V7 mp n lCT 

(ij)a ia ia 

+ V ^2n^n il + Y ^2 n i(T n j<7 r, (6) 
i {ij)<y<y' 

where (ij) stands for nearest neighbor pairs, and the 
notation is otherwise standard. Such a model, treated 
within a self-consistent BdG theory, has been used pre- 
viously to study vortices [p3|j2^ 1 and impurities |p^| , po| , p5| 
in d-wave superconductors. The impurities are modeled 
by V^ lmp = V lrop ^> \t\ at randomly chosen sites with den- 
sity rii and T/ Imp = elsewhere. We solve this Hamilto- 
nian within the standard mean-field theory as described 
in Ref. ^0|. All physical quantities of interest can be 
derived from the quasiparticle amplitudes _[u„(r), v n (jc)], 
which satisfy a system of BdG equations |2(| 
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£it n (r;) = -* y^u n (rj + 5) + (y/ mp - /i)un(ri), 

+ ^A 5 (r i )u„(r i + 5), (8) 

subject to the constraints of self-consistency 
A (r) = V J2un{r)v*(r)tanh(E n /2k B T), 

n 

A 5 (r) = ^ J] t u »( r + + w «( r X( r + *)] 

n 

x tanh(.E„/2fc B T). (9) 

Here Ao and As are the on-site and nearest-neighbor 
pairing amplitudes respectively with S = ±i, ±y and the 
summation is over positive eigenvalues E n only. For a 
finite L x L system we solve the BdG equations (0) by 
exact numerical diagonalization using a suitable guess 
for the initial gap functions. We then compute new gap 
functions from Eq. (^) and iterate this process until self- 
consistency is established |2^^] . 

The superfluid density p s (T) is evaluated from the 
standard linear response formula appropriate for lattice 
models (l§|27| 

Ps (T)/4 = (-K x ) - K xx {q x = 0, q y -> 0, w = 0), (10) 

where (—K x ) is the average kinetic energy along the in- 
direction, and A xx (q, u>) is a diagonal element of the 
current-current correlation function. Both (—K x ) and 
A xx (q, lu) can be written in terms of (u n ,v n ) and E n , 
resulting in lengthy expressions which we omit here for 
brevity. 

Determination of T c is more complicated since the 
number of iterations needed to self-consistently solve Eqs. 
(@iH) to a required accuracy increases dramatically near 
the transition (and presumably diverges at T = T c ). We 
use two different methods to overcome this difficulty: (i) 
we improve our initial guess for A(r) at each new temper- 
ature by extrapolating from previous temperature points, 
and (ii) we locate T c iteratively by following the devel- 
opment of the gap initialized to infinitesimal value. 

Figure || summarizes the results of our numerical cal- 
culation, for parameters resulting in £o ss 4ao (we use 
V = -Vi = 1.13*, (i = -0.36* and V imp = 100* for 
a system size L = 22) Panel (a) shows the normalized 
zero-temperature superfluid density as a function of rij 
obtained from a fully self-consistent solution of Eqs.(0J|) 
and from a solution with an enforced uniform order pa- 
rameter. For comparison a conventional analytic solution 
within the *-matrix formalism is also shown (we use Eq. 
(15) of Ref. §). As one would expect the uniform so- 
lution is very close to the analytical one for all values 
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FIG. 2. (a) Normalized zero-temperature superfluid den- 
sity as a function of impurity concentration: fully 
self-consistent solution of BdG equations (solid symbols), so- 
lution with uniform order parameter (open symbols) and an- 
alytic i-matrix solution with unitary scatterers from Ref. Q 
(dash-dotted line), (b) Normalized critical temperature and 
average gap. Dash-dotted line is numerical solution of the AG 
equation (Jl|). 

of rii. The self-consistent solution agrees well with the 
two, except for large m. Our results for T c [panel(b)] 
are much more interesting. While T c computed for the 
uniform solution tracks the analytic solution of the AG 
equation (|l|), the true critical temperature obtained us- 
ing a self-consistent solution is much higher. In fact at 
rii ~ 0.08 where T^ G — 0, the true critical temperature 
is still more than 70% of T c q. Also note that the aver- 
age zero-temperature gap does not scale with T c , as one 
would expect from BCS theory. 

We have carried out calculations for several differ- 
ent sets of parameters, all showing similar behavior (see 
Fig. ||). As the coherence length grows, the discrepancy 
between the true critical temperature and T^ G dimin- 
ishes, as expected from the analysis presented above. 
For £o/ao ~ 4 as in YBCO, the discrepancy is substan- 
tial, and can easily account for the experimentally ob- 
served robustness of T c . The critical impurity density 
for this case is n\ « 0.17. Since Zn substitutes primar- 
ily for the planar copper sites in YBCO [Q and only 
2/3 of all Cu atoms reside in the planes, our result im- 
plies a bulk critical density of n^ ulk w 0.10. This is in 
reasonable agreement with the experimentally observed 
n£ ulk = 0.08 - 0.10 for Zn doped YBCO ||. By con- 
trast conventional models tend to underestimate n£ ulk by 
a factor of 2 or more jl5| . The inset to Fig. || shows the 
dependence of the initial rate of change of T c with m, 
r]i = T^ ) 1 (dT c /dni) ni= o, on £o- This quantity also devi- 
ates significantly from the linear ru ~ (Co/^o) behavior 
expected from the AG prediction (ML showing instead a 
quadratic behavior consistent with (]J). 
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FIG. 3. Normalized critical temperature versus normalized 
zero-temperature supernuid density as computed from a fully 
self-consistent solution of BdG equations for systems with dif- 
ferent coherence lengths. The error bars reflect the statistical 
scatter of data for 6 different impurity configurations. Pa- 
rameters used are L = 22, F imp = -lOOt, fj, = -0.36t, and 
Vo = -Vi = (0.80, 1.05, 1.40)t for £ = (9.9, 4.7, 2.5)a re- 
spectively. Dash-dotted line is a conventional AG solution 
from Ref. Inset: the rate of change r/i versus £o/ao 

(open squares), and the expected behavior from AG theory 
r/i = 2.18(£o/ao) (dashed line). 



In closing we comment on the appropriateness of us- 
ing BdG theory to calculate T c in short coherence length 
materials. Like all conventional theories of superconduc- 
tivity, BdG theory is a mean field theory and hence is 
incapable of modeling the important effects of phase fluc- 
tuations on the ordering temperature in short coherence 
length, highly anisotropic materials such as the high-T c 
oxides. On the other hand, unlike many more specialized 
mean field theories, including AG theory, BdG theory is 
well-suited to calculating the effects of static spatial dis- 
order on the order parameter amplitude and on the mean 
field T c . Since the static random potential does not cou- 
ple to the phase, the effects of spatial randomness and 
thermal phase fluctuations are largely decoupled. To a 
good approximation, at least for the dilute case, the ef- 
fects of impurities can be calculated within mean field 
theory, while the effects of phase fluctuations, in reduc- 
ing T c from its mean field value, require a proper theory 
of critical behavior. A complete theory for the combined 
effects is not presently known. However it is apparent, 
even at the mean field level, that such a theory for the 
effect of inhomogeneous pair-breaking on T c must take 
proper account of the resulting inhomogeneity of the gap 
function. It is our hope that the present calculations will 
provide the physical motivation for such a theory. 
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